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We studied the statics and dynamics of elastic manifolds in disordered me- 
dia with long-range correlated disorder using functional renormalization group 
' (FRG). We identified different universality classes and computed the critical 

£^ ' exponents and universal amplitudes describing geometric and velocity-force 

Q characteristics. In contrast to uncorrelated disorder, the statistical tilt sym- 

O ■ metry is broken resulting in a nontrivial response to a transverse tilting force. 

For instance, the vortex lattice in disordered superconductors shows a new 
glass phase whose properties interpolate between those of the Bragg and Bose 
glasses formed by pointlike and columnar disorder, respectively. Whereas there 
is no response in the Bose glass phase (transverse Meissner effect), the stan- 
dard linear response expected in the Bragg glass gets modified to a power 
law response in the presence of disorder correlations. We also studied the long 
distance properties of the O(N) spin system with random fields and random 
anisotropics correlated as l/x d ~ a . Using FRG we obtained the phase diagram 
£f*) • in (d, <t, A r )-parameter space and computed the corresponding critical expo- 

f"*"^ ' nents. We found that below the lower critical dimension 4 + a, there can exist 

QQ ' two different types of quasi-long-range-order with zero order-parameter but 

, infinite correlation length. 
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H 1. Elastic manifolds in disordered media 

03 



Elastic objects in disordered media are a fruitful concept to study diverse 
physical systems such as domain walls (DW) in ferromagnets, charge den- 
sity waves (CDW) in solids, vortices in type-II superconductors. 1 In all 
these systems the interplay between elastic forces which tend to keep the 
system ordered, i.e., flat or periodic, and quenched disorder, which promotes 
deformations of the local structure, forms a complicated energy landscape 
with numerous metastable states. Most studies of elastic objects in dis- 
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ordered media restricted to uncorrclatcd pointlike disorder. Real systems, 
however, often contain extended defects in the form of linear dislocations, 
planar grain boundaries, three-dimensional cavities, etc. We studied the 
statics and dynamics of elastic manifolds in the presence of (power-law) 
long-range (LR) correlated disorder. 2 The power-law correlation of defects 
in d-dimensional space with exponent a = d — Sd can be ascribed to ran- 
domly distributed extended defects of internal dimension Sd with random 
orientation. 

The configuration of elastic object can be parametrized by displacement 
field u Xl where x denotes the d dimensional internal coordinate of the elastic 
object. The Hamiltonian is given by 



c is the elasticity, V(x, u) a random potential with zero mean and variance 



V{x,u)V{x',u') = Rt{u - u')5 d {x - x') + R 2 (u - u')g{x - x'), (2) 



where g(x) ~ x~ a . For interfaces one has to distinguish two universality 
classes: random bond (RB) disorder described by a short-range function 
R(u) and random field (RF) disorder corresponding to a function which 
behaves as R(u) ~ |u| at large u. Random periodic (RP) universality class 
corresponding to a periodic function R(u) describes systems such as CDW 
or vortices in d — 1 + 1 dimensions. Disorder makes the interfaces rough 
with displacements growing with the distance x as C(x) ~ a; 2 '', where C is 
the roughness exponent. Elastic periodic structures lose their strict trans- 
lational order and may exhibit a slow logarithmic growth of displacements, 
C(x) = Ad In | x | . The driven dynamics at zero temperature is described by 
the over-damped equation of motion 



Here r] is the friction coefficient, / the applied force and F = —d u V(x, u) 
the pinning force with zero mean and correlator 



F(x, u)F(x',u') =A 1 (u- u')S d (x - x') + A 2 (u - u')g(x - x'), (4) 



where Ai(u) = —R'l(u). The system undergoes the so-called depinning 
transition at the critical force / c , which separates sliding and pinned states. 
Upon approaching the depinning transition from the sliding state / — > /+ 
the center-of-mass velocity v = L~ d j dtu xt vanishes as a power law v ~ 
(/ ~ /c) and the time and length scales are related by t ~ x z . 




(1) 



ildtu x t 



c\7 2 u xt + F{x,u xt 



) + /■ 



(3) 
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Let us firstly discuss the case of short range (SR) correlated disorder 
(A 2 =0). The problem turns out to be notably difficult due to the so-called 
dimensional reduction (DR). However, simple Imry-Ma arguments show 
that DR is wrong. The metastability renders the zero temperature pertur- 
bation theory useless breaking it down on scales larger than the so-called 
Larkin length. To overcome these difficulties one can apply functional renor- 
malization group (FRG). Scaling analysis gives that the large-scale proper- 
ties of a d-dimensional elastic system are governed by uncorrelated disorder 
in d < d uc — 4. The peculiarity of the problem is that there is an infinite set 
of relevant operators which can be parametrized by a function. The latter 
is nothing but the disorder correlator. We generalized the FRG approach to 
systems with LR correlated disorder. Using a double expansion in e = 4 — d 
and 8 = 4 — a we derived the flow equations to one-loop order: 

deA^u) = (e- 2()A 1 (u) + (uA[(u) - ^[Ai(«) + A 2 (u)] 2 + AA'{{u), 

d e A 2 (u) = (5- 2C)A 2 («) + (uA' 2 (u) + AA' 2 '{u), (5) 

where A = [Ai(0) + A2(0)]. The scaling behavior of the system is controlled 
by a stable fixed point [A*(u), A 2 (u)] of flow equations (5). Note that for 
uncorrelated disorder the elasticity remains uncorrected to all orders due 
to the statistical tilt symmetry (STS). LR disorder destroys the STS and 
allows for the rcnormalization of elasticity introducing a new exponent ip. 
The SR part of disorder correlator Ai (it) becomes a nonanalytic function 
beyond the Larkin scale, while the LR part A 2 (u) remains analytic along 
the flow. The appearance of non-analyticity in the form of a cusp at the 
origin is related to metastability, and nicely accounts for the generation of 
threshold force f c at the depinning transition. The actual values of critical 
exponents are fixed by the disorder correlator at the stable FP, for instance, 
tjj = \(5 - e)A*"(0) and z = 2 - A*"(0) - A*"(0). Let us summarize the 
results obtained for different universality classes. 

RB disorder. The crossover from the SR FP (A 2 = 0) to the LR FP 
(A 2 ^ 0) takes place for 5 > 1.041e. At the LR RB FP there is an exact 
relation between exponents: £lr = (S + 2^>)/5, where ip — 0(e 2 ,5 2 ,e8) to 
one-loop order. We also computed the amplitude of height-height correla- 
tion function which turns out to be universal up to the strength of disorder. 

RF disorder. The large-scale properties of the system is controlled by the 
LR FP for S > e. The roughness exponent is given by the exact relation 
Clr = (S + 2tp)/3. Surprisingly, the function Al(u) satisfies du A\(u) = 
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0, characteristic for RB type correlations along the u direction. In other 
words, the LR RF FP is in fact of mixed type: RB for the SR part and 
RF for the LR part of the disorder correlator. The RF FP describes the 
depinning transition of interfaces in the presence of LR correlated disorder. 
The corresponding dynamic exponents read: z = 2 — e/3 + (5/9 and (3 = 
1 — e/6 + 5/18. It is remarkable that for 8 > 3t~ the exponent (3 is larger 
than 1, and z larger than 2. 

Random periodic. For 5 > n 2 s/9 the periodic system is controlled by 
the LR RP FP. The system exhibits a slow growth of the displacements: 
(u x — u ) 2 — In \x\. As an example, we consider vortices in the presence 
of three types of disorder: uncorrelated, LR and columnar. The correlation 
of disorder violates the STS resulting in a highly nonlinear response to tilt. 
In the presence of columnar disorder vortices exhibit a transverse Meiss- 
ner effect: disorder generates the critical field h c below which there is no 
response to tilt and above which the tilt angle behaves as d ~ (h — h c )^ 
with the universal exponent (f> < 1. The RP case describes a weak Bose 
glass which is expected in type-II superconductors with columnar disorder 
at small temperatures and at high vortex density which exceeds the density 
of columnar pins. The weak Bose glass is pinned collectively and shares fea- 
tures of the Bragg glass, such as a power-law decay of the translation order, 
and features of the strong Bose glass such as a transverse Meissner effect. 
For isotropically LR correlated disorder the linear tilt modulus vanishes 
at small fields leading to a power-law response # ~ with (f> > 1. The 
response of the system with LR correlated disorder interpolates between 
responses of systems with uncorrelated and columnar disorder. We argued 
that in the presence of LR correlated disorder vortices can form a strong 
Bragg glass which exhibits Bragg peaks and vanishing linear tilt modulus 
without transverse Meissner effect. 3 

2. LR random field and random anisotropy O(N) models 

The large-scale behavior of the O(N) symmetric spin system at low tem- 
peratures can be described by the nonlinear a model with the Hamiltonian 



where s(x) is the iV-component classical spin with a fixed-length constraint 
s 2 = 1. V(x, s) is the random disorder potential, which can be expanded in 



spin variables: V{x, s) = - J2™=i £n...i„ h^> l Jx)s il (x)...s i ^x). The cor- 
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responding coefficients have simple physical interpretation: h^p is a random 
field, fiff is a random second-rank anisotropy, and are general fith ten- 
sor anisotropies. We studied system (6) with LR correlated disorder given by 

cumulants/^.^x)^..^')- 5^5 iljl ...6 i ^[r^8{x-^)+r^g{x-x')] 
with g(x) ~ |x| <J_d . The corresponding replicated Hamiltonian reads 4 



H n = J d d .r 



2T 

a.b 



(7) 



- ^ E/ d<V -^ - x') R 2 (s a (x) ■ s b (x')) 

a.b 

where Ri{z) — X^ 7 ^ 2 ^ are arbitrary for RF and even for RA. Power 
counting suggests that d\ c = 4 + a is the lower critical dimension for both 
models. At criticality or in the quasi-long-range-ordered (QLRO) phase the 
connected and disconnected two-point function scale with different expo- 
nents i] and fj. The FRG equations to first order in e = 4 — d and o are 
given by 

d t Ri{<t>) = -eR 1 {4>) + \[R'l{4>) + R'^)] 2 — (N — 2)^2AR 1 {4>) 

-AR'{(cj)) + AR[ (</>) cot - — L- + R' 2 (<j>)Y 

2 sin <p 

d t R 2 {(j>) = -{e-<y)R 2 {4>)-[{N -2)[2R 2 {4>) + R' 2 {4>) cot 4>] +R'^)}a. 

Here ^4 = R"(0) + R 2 (0) and z = cos<f>. An attractive FP of the flow 
equations describes a QLRO phase, while a singly unstable FP describes 
the critical behavior. The critical exponents are determined by the FRG 
flow in the vicinity of the FP and to one-loop order are given by rj = —A* 
and fj = — e — (N — 1)A* . A singly unstable FP has only one positive 
eigenvalue X\ which determines the third independent exponent v = 1/Ai 
characterizing the divergence of the correlation length at the transition. For 
the RF model the critical exponents are given by 

£ - a _ 2s - (N - 1W 

The stability regions of different FPs and the phase diagrams for the RF 
and RA models are shown in Figure 1. For the LR RF model we found that 
the truncated RG developed in work 5 for the RF model gives the correct 
one- loop values of exponents rj and fj above the d\ c , but not the phase 
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Fig. 1. (a) Stability regions of various FPs for RF O(N) model corresponding to dif- 
ferent patterns of the critical behavior above the lower critical dimension, e > cr. In- 
sert: Schematic phase diagram on the (e, cr) plane for a particular value of N 6 (3, 18). 

(b) Phase diagram of the RF model below the lower critical dimension, d < 4 + a. 

(c) Phase diagram of the RA model below the lower critical dimension, d < 4 + cr. 

diagram and the critical exponent v except for the region controlled by the 
weakly nonanalytic LR TT FP. Thus, although the truncated RG overcomes 
the dimensional reduction it fails to reproduce all properties which can be 
obtained using FRG. We found a new LR QLRO phase existing in the 
LR RF model below the d\ c for N < 3 and determined the regions of its 
stability in the (e, cr, N) parameter space. We obtained that the weak LR 
correlated disorder does not change the critical behavior of the RA model 
above d\ c for N > N c — 9.4412, but can create a new LR QLRO phase 
below d\ c . The existence of two QLRO phases in the LR RA system may 
be relevant to two different states of 3 He-A in aerogel observed recently in 
NMR experiments. 6 
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